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als exactly coincide with orthogonal polynomials appearing in Lanczos tridiagonalization 
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1 Introduction 



Recently Bender and Dunne showed that the solution 'iIj{x, E) of the Schrodinger equation 
Hip{x, E) = Eip{x, E) for the one-dimensional quasi-exactly solvable model with hamiltonian 
H = —d'^jdx^ + c(c — + (4M + 2c + 3)x^ + is the generating function for the set of 

polynomials Pn{E) in the energy variable E. These polynomials, defined as the coefficients 
of the expansion of function exp(x^/4)'0(x, -E) in powers of x^, satisfy the three-term 
recursion relations and therefore form an orthogonal set with respect to some weight function 
ijj{E). Bender and Dunne demonstrated that a) for non-negative integer values of M for which 
the model K is quasi-exactly solvable, the norms of polynomials PniE) vanish for n > M + 1; 
b) the zeros of the critical polynomial PM+iiE) coincide with the exactly calculable energy 
levels in model H; c) the polynomials PniE) of degrees higher than M + 1 factor into a 
product of two polynomials, one of which is Pm+i{E): PM+i+n{E) = PM+i{E)Qn{E); d) 
the factor-polynomials Qn{E) also form an orthogonal set. 

These and other remarkable properties of polynomials Pn{E) discussed in paper |Q| sug- 
gest that they can be considered as an alternative mathematical language for describing the 
phenomenon of quasi-exact solvability. This phenomenon has been discovered several years 
ago (see e.g. the original papers ^, |^, |5| and books 1^, 0]). In their paper, Bender and 
Dunne did not present any proof of the existence of polynomials Pn{E) with similar properties 
for other quasi-exactly solvable models. Many important questions, such as the problem of 
construction of some analogues of Bender - Dunne polynomials for multi-dimensional quasi- 
exactly solvable models, the problem of the computation of weight functions associated with 
polynomials Pn{E) and the corresponding factor-polynomials Qn{E), also remained to be 
open. 

The aim of the present paper is threefold. First of all, we intend to propose a model inde- 
pendent way of introducing Bender - Dunne polynomials for quasi-exactly solvable systems 
of quantum mechanics. We show that the polynomials PniE) (which we call the Bender - 
Dunne polynomials of genus one) and the corresponding factor-polynomials Qn{E) (which 
we call the Bender - Dunne polynomials of genus two) naturally appear as the coefficients of 
the expansion of solutions %l){x,E) of Schrodinger equation in the bases in which the hamil- 
tonian H has a tridiagonal form Q Second, we present a simple way for building the weight 
functions uj{E) and p{E) associated with polynomials Pn{E) and ] Qn{E) and show that 
these polynomials are orthogonal polynomials of a discrete variable E^. In the case of PniE) 
polynomials, the variable E^ takes its values in a finite set of the exactly calculable energy 
levels in model H, and, in the case of Qn{E) polynomials, it runs an infinite set of all the 
remaining (exactly non-calculable) levels. Third, we demonstrate that the Bender - Dunne 
polynomials exactly coincide with orthogonal polynomials in an operator variable H arising 
in the so-called Lanczos tridiagonalization procedure [§, |^, 10]. Since this procedure is ap- 



plicable to any quasi-exactly solvable models (both one-dimensional and multi-dimensional), 
this enables one to make a conclusion of the universality of Bender - Dunne construction. 



''Note that the hamiltonians of all quasi-exactly solvable models are explicitly tridiagonalizable (for details 
see section 7-9) 
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2 Bender — Dunne polynomials of genus one 



Consider a hamiltonian H acting in Hilbert space W in which it has an infinite and discrete 
spectrum Ek-, fc = 0, 1, . . . , oo. Let re = 0, 1, . . . , oo be a certain (not neccessarily orthog- 
onal) basis in W in which the hamiltonian takes an explicit tridiagonal form. This means 
that 

H(j)n= An(t)n-l+ BnCkn + CnCkn+l-, n = 0, 1,...,00 (2.1) 

where and Cn are some algebraically computable coefficients. Note that values of 

these coefficients depened also on the normalization of basis vectors (pn- 
Assume that there exists a certain non-negative integer M for which 

Cm = 0, Cn / 0, for re / M. (2.2) 

Then, the model with hamiltonian H becomes quasi-exactly solvable. Indeed, from formulas 
( |2.1| ) and ( |2.2| ) it follows that a linear hull of the basis vectors (/>o, . . . , (I^M forms a (M -|- 1)- 
dimensional invariant subspace Wm C W for the hamiltonian H. For this reason, the spectral 
problem for H m W breaks up into two independent spectral problems, one of which is 
formulated in the (M -|- l)-dimensional space Wm- This enables one to determine at least 
M + 1 eigenvalues Eq, Ei, . . . , Em of hamiltonian if in a purely algebraic way. The remaining 
eigenvalues Em+i, Em+2, ■ ■ ■ cannot be determined algebraically and therefore the model 
under consideration is quasi-exactly solvable (see ref.|^). Hereafter we shall call the levels 
£^0) • • • ) Em exactly calculable and the levels Em+i, Em+2, ■ ■ ■ - exactly non-calculable. 
Let us now consider the Schrodinger equation 

H^p{E) = E^P{E) (2.3) 

in which E is an arbitrary parameter. We look for its solution in the form of the formal 
expansion 

oo 

^{E) = ^ ^nPn{E) (2.4) 
n=0 

in which Pn{E) denote certain functions of E. 

Note that the expansion (^]J) diverges for E ^ E^ and converges for E = E^. If the 
model is quasi-exactly solvable, this is, if the condition ( p.2| ) is satisfied, then the series (2.4) 
becomes finite for all E = Ef^, k = 0, . . . ,M. This means that for these values of E the 
values of functions Pn{E) should vanish for re > M + 1: 

Pn{Ek) = Q, re>M + l, A: = 0,1,...,M. (2.5) 

In order to determine the form of functions Pn{E) it is sufficient to substitute the ex- 
pansion ( |2.4D into equation ( |2.3| ) and use the tridiagonality condition ( p.l[ ). This gives the 
recurrence relations for Pn{E), 

EPniE)=An+lPn+liE)+BnPniE)+Cn^lPn-liE), re = 0, 1 , . . . , OO, (2.6) 

which can be supplemented by the initial conditions 

P_i(E) = 0, Po{E) = l. (2.7) 
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From formulas (|2.6| ) and (^3) it immediately follows that functions Pn{E) are polynomials 
of degree n. From the general theory of orthogonal polynomials (see e.g. we know that 
Pn{E) are orthogonal polynomials. This means that there exists a certain weight function 
uj{E), which can be normalized as 

J ujiE)dE = 1, (2.8) 

for which 

J Pn{E)P^{E)oo{E)dE = PnPmSnm, (2.9) 

where pn denote the norms of polynomials Pn{E). These norms can be found from the 
recurrence relations ( p.6| ) after multiplying them by E'^~^lj(E) and taking integral over E. 
The result is 

1 

(2.10) 

If the model is quasi-exactly solvable, then from the condition ( |2.2D and formula ( |2.10| ) it 
follows that 

Pn = 0, n>M + l, (2.11) 

i.e. the norms of all polynomials Pn{E) with n > M + 1 vanish. Hereafter, we shall call the 
polynomials Pn{E) the Bender - Dunne polynomials of genus one. 



Pn 



m=l "^"^ 



3 Factorization property 

The Pn{E) polynomials of degrees n> M + 1 have a remarkable factorization property first 
noted by Bender and Dunne ||l|. Indeed, from formula (|2.5| ) it follows that the M + 1 exactly 
calculable eigenvalues Eq, Ei, . . . , Em of hamiltonian H are the roots of these polynomials. 
This enables one to write 

M 

PM+liE) ^ YliE - En) (3.1) 

n=0 

and 



P 



M+l+n 



{E) = PM+l{E)Qn{E), 



(3.2) 



where Qn{E) are certain polynomials of degrees n = 0, 1, . . . , oo. 

The first formula (3^) gives a practical way of determining the exactly calculable eigen- 
values of the hamiltonian H from the equation 



Pm+i{E) = 0. (3.3) 

The second formula ( |3.2[) expresses an important factorizability property of the polynomials 
Pn{E) with n > M + 1. The factor-polynomials Qn{E) we shall call the Bender - Dunne 
polynomials of genus two. 
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4 Bender — Dunne polynomials of genus two 



Substituting (3.2) into ( |2.6D , one can obtain the recurrence relations immediately for the 
factor-polynomials Qn{E): 

EQn{E) = AM+2+nQn+l{E) + BM+l+nQn{E) + CM+nQn-l{E), n = 0, 1, . . . , OO. (4.1) 

Since Cm = 0, it is sufficient to supplement the relations (4.1) by the initial conditions 

Qo{E) = 1. (4.2) 

We see that Qn{E) are again the orthogonal polynomials. Denoting the corresponding weight 
function by p{E) and normalizing it as 



j p{E)dE = 1, 



we can write 



QniE)Qm{E)p{E)dE = 



The norms of polynomials Qn{E) computed from the recurrence relations (4.1) are 



n ^ 
■j-j- L^M 



m=l 



(4.3) 



(4.4) 



(4.5) 



According to (p.2|), the norms of Qn{E) polynomials do not vanish. 

5 Weight functions 

In this section we compute the weight functions u){E) and p{E) for Bender - Dunne orthogonal 
polynomials of genus one and two. 

Weight function for Pn{E) polynomials. From formulas ( p.TD - ( |2.9[ ) it follows that 



j Pn{E)uj{E)dE = 6nO. (5.1) 

Multiplying the relation ( p.4[ ) by uj{E), taking the integral over E and using ( ^.1[) , we obtain 

ilj{E)uj{E)dE = (j)o. (5.2) 



The basis element (po can obviously be represented as a linear combination of eigenvectors of 
hamiltonian H corresponding to exactly calculable energy levels: 



M 



bo = ^iVkipiEk). 

k=0 



(5.3) 



Substituting ( ^.3[ ) into (5^), we obtain 

. M 
J il,{E')uj{E')dE' = Y,uJkil>(Ek). 



(5.4) 



fc=0 
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Let us now act on both hand sides of the relation ( p.4|) by the operator function 5{H — E). 
Using (|2.3D , we obtain 

M 

^{E)u{E) = Uk5{E - EMEk). (5.5) 

fc=0 



Dividing both hand sides of (5^) by ■0(-^) S'lid taking into account the projection property 
of delta- function, we obtain the final expression for the weight function uj{E): 

M 

oj{E) = Y,oJk5{E-Ek). (5.6) 

fc=0 

Thus, we have arrived at interesting result: the Bender - Dunne polynomials of genus one 
are orthogonal polynomials in a discrete variable Ek, A; = 0, 1, . . . , M which takes its values 
in the finite set of the exactly calculable eigenvalues of hamiltonian H. The orthogonality 
property for these polynomials can now be written in the form 

M 

(5.7) 

k=0 



The coefficients lo^ can be computed algebraically. Indeed, substituting into ( ^.31) the expan- 
sions (2.4), we obtain the system of algebraic equations for numbers uok'- 

AI 

J2 PniEk)iOk = 6no, n = 0,l,...,M. (5.8) 

fc=0 

The numerical analysis of this system for various models shows that the numbers cok, playing 
the role of a discrete weight function are always positive (see e.g. section 6). It would be an 
interesting task to proof this in general case. 



Weight function for Qn{E) polynomials. Using formulas (|4.2|) - (|4.4| ) we obtain 

J Qn{E)p{E)dE = 6no. (5.9) 
Using the factorization property (|3.2[), let us represent the expresion (|2.4D in the form 



A/ 



V^(^) = J2 ^nPn{E) + Pm+i{E) <pM+l+nQniE). (5.10) 
n=0 n=0 

Multiplying ( ^.10| ) by p{E)Pf^j^^^{E), taking the integral over E and using (|5.9D, we obtain 

J Fm+i[E') J Fm+i[E') 

Expressing vectors (/>„, n = 0, 1, . . . ,M and (t)M+i via the orthonormalized eigenvectors of 
hamiltonian H, we can rewrite ( ^.ll| ) as 



/ J^^'l,M E')dE' = J2 f^ki^iEk) + E i^ki^iEk), (5.12) 
J Pm+i{E) , „ 
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where fik and are certain, in general, algebraically non-computable coefficients. Note that 

i^k = {<pM+i,HEk)), k = M + l,...,oo, (5.13) 

where ( , ) denotes the scalar product in Hilbert space W. Let us now act on both hand 
sides of the relation ( |5.12D by the operator function 6{H — E). Using (|2.3| ), we obtain 

' -p{E) = J2fik6{E-EMEk)+ E ^kS{E-Ek)i:{Ek). (5.14) 



PM+iiE)' 



k=0 k=M+l 
j-1 



Dividing both hand sides of ( 5.14 ) by 'ip{E)P^j^-^{E) and taking into account the projection 
property of delta- function and formula ( |2.5| ), we obtain the final expression for the weight 
function p{E): 

oo 

piE)= PkSiE-Ek), (5.15) 

k=M+l 

where 

Pk = PM+i{Ek)i^k = PM+i{Ek){4>M+i,'4^{Ek)), k = M + 1, . . . ,oo. (5.16) 

Thus, we see that the Bender - Dunne polynomials of genus two are the orthogonal polyno- 
mials in a discrete variable E/., k = M + 1, . . . ,oo whose values are nothing else than the 
exactly non-calculable eigenvalues of hamiltonian H. The orthogonality property for these 
polynomials can now be written in the form 

oo 

PkQn{Ek)Qm{Ek) = (5.17) 

k=M+l 

Note that the numbers pk cannot be found by means of algebraic methods. 



6 An example 

In this section we construct the discrete weight function ojk for Bender - Dunne polynomials 
of genus one associated with the simplest quasi-exactly solvable models with polynomial 
potentials. The hamiltonians of these models, which were first discussed in ref. Q, have the 
form 



H 



+ + 26x^ + [5^ - 4M - 3]a;2, 



(6.1) 



and for any given non-negative integer M admit M + 1 explicit solutions corresponding to 
energy levels Eq,Ei,. . . ,Em describing the first M + 1 even states with ordinal numbers 
0,2,..., 2M. It is easily seen that in the basis 



2\n 



{2n)l 



■ exp 



X 

'T 



bx^ 
2 



(6.2) 



the hamiltonian H takes an explicit tridiagonal form. The corresponding coefficients An, Bn 
and C„ are 



An = l, Bn = b{4n + 1), Cn = 8a{n + l){2n + l){n- M). 



(6.3) 
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The general formulas of sections 2-5 allow one to reconstruct the polynomials Pn{E) and 
compute the corresponding discrete weight functions. Below we consider two examples of 
such calculations. 

Example 1. Let us take 6 = and M = 8. Then Bender - Dunne polynomials of genus 
one reconstructed from recurrence relations read 

P^[E) = 1, 

Pi{E) = E, 

P^(E) = -64 + 

Ps{E) = -mE + E^, 

Pi{E) = 46080 - 1120 E"^ + E^, 

P5 {E) = 494080 E - 2240 E^ + E^ , 

Pq{E) = -66355200 + 2106880 E'^ - 3680 E'^ + S^ 

Pr^[E) = -848977920 E + 5655040 E^ - 5264 E^ + E\ 

Ps{E) = 96613171200 -3916595200 + 11013120 -6720 +^^, 

Pg{E) = 911631974400 E - 9345433600 E^ + 16066560 E^ - 7680 E'^ + E^ . (6.4) 

The energy levels coinciding with the roots of the last polynomial Pq{E) are 

Eq = -68.1568, 

El = -46.5609, 

E2 = -27.2997, 

^3 = -11.021, 

Ea = 0, 

^5 = 11.021, 

Eq = 27.2997, 

Ej = 46.5609, 

^8 = 68.1568, (6.5) 

and the components of the corresponding weight function oo^ computed by means of formula 
( |5^ ) are 

loq = 0.0000117487, 

= 0.00058638, 
UJ2 = 0.013045, 
^3 = 0.172498, 
UJA = 0.627717, 

= 0.172498, 
We = 0.013045, 

= 0.00058638, 
a;8 = 0.0000117487. 

We see that all these numbers are positive. Note also that these numbers are disctributed 
symmetrically with respect to the reflection /c — > M + 1 — fc, while the energy levels are 
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distributed anti-symmetrically. This is a trivial consequence of the self duality property of 



the model ( |6.1| ) with 6 = noted in ref. |12|. 



Example 2. Let us now take b = 3 and M = 8. Then Bender - Dunne polynomials of 
genus one reconstructed from recurrence relations read 

Po{E) = 1, 
Pi{E) = -3 + E, 
P^[E) = -19-18^ + ^^ 
Pz{E) = 1521 + 131 45 £^^ + ^3, 
P^{E) = -45639 + 9372 £ + 1166^2- 84 E^^ + E^^, 
P^{E) = 624069 - 670331 £ + 306 + 4330 135 ^;^ + ^;^ 
Pq{E) = 26403813 + 29359242 £- 2368649 £^2- 151524 £^3 + 11395 + 
-198£^ + £^ 

Prji^E) = -2968811271 - 1113735033 E + 206523213 £^ + 2136931 E^ + 

-792309 E^ + 24661 E^ - 273 E^ + E\ 
Ps{E) = 219842628849 + 51179080248 £ - 15632501620 £^ + 241229160 £^ + 

+54476694 £^ - 2649528 £^ + 46956 £^ - 360 £'^ + £^, 
Pg{E) = -18914361435891 - 3777700684023 £ + 1400534456148 E^ + 

-41565642220 £3 - 4391346906 £^ + 293105406 £^ - 7036092 £6 + 

+81636£^ -459£'^ + £^. (6.6) 

For the energy levels we have 



£o 


= -14.4717 


El 


= -2.52391 


£2 


= 6.59852, 


£3 


= 21.606, 


£4 


= 40.6664, 


E5 


= 62.6983, 


Ee 


= 87.2765, 


£7 


= 114.122, 


£« 


= 143.028, 



(6.7) 



and the components of the weight function are 



Wo = 0.026295, 

coi = 0.475485, 

u;2 = 0.423177, 

W3 = 0.067025, 

a;4 = 0.00741457, 

W5 = 0.000573463, 

loq = 0.0000291886, 
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= 8.80537 lo- 
ws = 1.1954 10"^ 



(6.8) 



We see that the discrete weight function is again positive. At the same time, we have no 
symmetry of its components and no anti-symmetry of the energy levels, because for b ^ 
the model (|0]) is not self-dual |12]. 



7 One-dimensional quasi-exactly solvable models 



The hamiltonians H = —d'^jdx^ -|- V{x) of one-dimensional quasi-exactly solvable models 
introduced in paper [Q] are distinguished by the fact that for them it is possible to find two 
functions A(x) and p{x) for which the ansatz 



reduces the corresponding Schrodinger equation H'\\){x^E) = Eip(xjE) to the form 



o2 f) 

i?3(A)^ + i?2(A)^+i?i(A) 



f{X,E) = EfiX,E) 



(7.1) 



(7.2) 



with some polynomials i?i(A), R2{X) and -R3(A) of degrees mi = 1, 1712 
respectively. The phenomenon of quasi-exact solvability appears when 

M(M- 1) lim X~^R3{\)+M lim A"2i?2(A) + lim A"^i?i(A) : 

A^oo A^oo A^oo 



2 and ms < 3, 



0, 



(7.3) 



where M is a certain non-negative integer. It is easily seen that, for any M, there exist M + 1 
values of E = Ek for which the equation ( fT^ ) has solutions f{X,E) = f{X,Ek) represented 
by certain polynomials of degree M. From ( |7.3| ) it immediately follows that in the basis 

(7.4) 



(pnix) = [X{x) - r]'^p{x), n = 0, 1, . . . ,cx) 



in which r denotes any of the roots of the polynomial i?3(A), the hamiltonian of the model 
(irrespective of whether the condition of quasi-exact solvability is satisfied or not) takes an 



explicit tridiagonal form (2.1) with trivially computable coefficients A^, Bn and C„. Indeed, 
introducing the new variable t = X — r and taking f{X,E) = g{t,E), we obtain instead of 

(ItT 



d 



d 



(03^ + b^t^ + C3t)-g^ + (aar + b2t + 02)7^ + (ait + 61) 



dt 



git,E) = Eg{t,E). (7.5) 



It is easily seen that the operator L standing in the left-hand of the equation ( [7.5[ ) is a linear 
combination of three operators L"*", L^, L~ acting on the monomials as L'^t'^ ~ t 



n+l 



L^t"^ ~ and L ^ t"^ ^. But this means that this L-operator is tridiagonal in the basis 



t", n = 0, . . . , 00. This, in turn, results in the tridiagonality of H in basis (7.4). According to 
general prescriptions given above, this means that any one-dimensional quasi-exactly solvable 
model from the list given in ref. |Q] trivially generates ms < 3 Bender - Dunne polynomials 
of genus one and two. An explicit construction of these polynomials and computation of the 
corresponding weight functions will be given in a separate publication. 
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8 Lanczos tridiagonalization procedure 



There are infinitely many ways of reducing a given hamiltonian if to a tridiagonal form. One 
of the simplest ways is based on the use of the so-called Lanczos tridiagonalization procedure 
(see e.g. refs.||8|, ^, |^]). In terms of the Lanczos procedure, the properties of Bender - Dunne 
polynomials become especially transparent and simple. 

Let -ff be a hamiltonian and </> € be an arbitrary vector of Hilbert space. Consider the 
sequence of vectors 

(j)n = Pn{H)^, n = 0,l,...,oo, (8.1) 

in which Pn{H) denote certain polynomials in the operator variable H. Let us choose these 
polynomials from the condition of orthonormalizability of vectors (pn 

{4>n,4'm) = Snm (8.2) 

with respect to the scalar product in the space W. If the hamiltonian is given and the scalar 
product is explicitly defined, then the orthogonalization is a purely algebraic procedure. The 
sequence 0„ can be considered as an orthogonal basis in Hilbert space. Let V'C-^fe) be the 
orthonormalized eigenvectors of the hamiltonian H, 

{il,{Ek),il^{Ei)) = 6u. (8.3) 
Expanding the vector (j) in eigenvectors iplE/^), 

oo 

= 5]cfcV(Sfc), (8.4) 

k=0 



and substituting this expansion into ( |8.1| ) we obtain 

oo 

cl>n = Y.''kPn{EMEk). (8.5) 

k=0 

The substitution of ( |8.5| ) into the orthonormalizability condition (^]^) gives the relation 

oo 

J2 clPniEk)Pm{Ek) = Snm (8.6) 
A:=0 

which shows that Pn{E) are orthogonal polynomials in a discrete variable Ek^ A; = 0, 1, . . . , oo, 
which takes its values in the set of eigenvalues of hamiltonian H (see ref . [lTo|| ) . The corre- 
sponding discrete weight function c| is obviously positive. Hereafter, we shall call Pn{H) 
the Lanczos polynomials. Because of the orthogonality, the polynomials Pn{H) obey the 
recurrence relation 

HPn{H)= AnPn^l{H)+BnPn{H)+CnPn+l{H), n = 0, 1, . . . , OO, (8.7) 

in which Bn and C„ are certain algebraically computable coefficients. Acting by the 



operator equality (8^) on (j) we obtain 



H(l)n = AnCt)n-l + Bn4>n + Cn4>n+1, n = 0, 1, . . . , OO, 
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which means that the operator H is tridiagonal in the basis (j)n, n = 0, 1, . . . , oo. Since H 
is a hermitian operator and the basis is orthogonal, the upper and lower diagonals of H 
should coincide. This gives us the conditions 

Cn = (8.9) 

which reduce the recurrence relations (|8.7D to the form 

HPn{H)=Cn-lPn-l{H)+BnPn{H)+CnPn+l{H), n = 0, 1, . . . , CX). (8.10) 

Let us now consider the equation (^]^) for the hamiltonian H. Looking for its solutions in 
the form (|2^), substituting (0) into (|8^ and using formula ( |8.9D we find that Bender - 



Dunne polynomials P„ associated with the expansion ( p.^ ) satisfy the recurrence relation 

EPn{E)=Cn-lPn-l{E)+BnPn{E)+CnPn+l{E), n = 0, 1, . . . , OO. (8.11) 



From the coincidence of the recurrence relations ( 8.10| ) and ( ^.11 ) it follows that the Bender 



- Dunne polynomials Pn{E) are nothing else that the Lanczos polynomials Pn{E). Hereafter, 
we shall not distinguish between the Lanczos and Bender - Dunne polynomials and for both 
of them use the same notation Pn{E). 

Up to now, we implicitly assumed that all the coefficients cj. in the expansion (|8.4| ) differ 
from zero. Assume now that this expansion has only a finite number of terms 

M 

cj)=J2ckHEk). (8.12) 

k=0 

Let us denote by Wm the linear span of functions ijj{Ek), A; = 0, 1, ... , M. The space Wm is 
a (M + l)-dimensional invariant subspace for the hamiltonian H. It is quite obvious that, in 
this case, the Lanczos procedure becomes finite and consists only of M + 1 essential steps. 
Indeed, assume that the Lanczos polynomials Pn{E) with < n < M are already known and 
we want to construct the next polynomials PM+i+n{H). According to general prescriptions 
of Lanczos theory, the vectors cjiM+i-^n should have the form (f)M+i+n = PM+i+n{H)4> and 
thus should belong to the space Wm- On the other hand, they should be orthogonal to all 
the linearly independent basis vectors (po, (pi, . . . , (pM of the space Wm, which means that 
<pM+i+n = PM+i+n{H)(j) = 0. This is possible only if PM+i+n{H) = ]\ILq{H - Ek)Qn{H), 
where Eq, Ei, . . . , Em are the eigenvalues of the hamiltonian H in the space Wm and Qn{H) 
are the arbitrary polynomials of degrees n = 0, 1, . . . , oo. The fact that the functions (j)M+i+n 
are non-normalizable, prevents one from determining the polynomials Qn{H) uniquely. 
The only way to construct the polynomials Qn{H) is to choose a new vector 

oo 

cPm+i = Pm+i{H) CkHEk) (8.13) 

k=M+l 

lying outside the space Wm and apply to it the Lanczos tridiagonalization procedure. This 
vector is, obviously, normalizable, and starting with it, we can construct the set of orthonor- 
malized vectors (pM+i+n = Qn{H)4'M+i belonging to an orthogonal complement of the space 
Wm- This procedure enables one to determine the polynomials Qn{H) uniquely. 

It is easily seen that the orthonormalizability condition for functions (pm n = 0, 1, . . . , M 
reads 

M 

J2 clPniEk)Pm{Ek) = Snm (8.14) 
k=Q 
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and is nothing else than the orthonormahzabiUty condition for Bender - Dunne polynomials 
of genus one. Note the positive definiteness of the corresponding weight function Wfc = c|. 
At the same time, the orthonormalizability condition for functions (l)M+i+n, = 0, 1, . . . , oo 
takes the form 

M 

^ clPlj+^iEk)QniEk)QmiEk) = 5nm (8.15) 
k=0 

and becomes that for Bender - Dunne polynomials of genus two. The corresponding weight 
function pk = c\Pljj_i{Ek) is also positively defined. 



9 Conclusion 

In conclusion note that the Lanczos tridiagonalization procedure can be applied to any her- 
mitean operators in Hilbert space. This means that H may be the hamiltonian of an ar- 
bitrary one or multi-dimensional model. However, if we want to construct the Bender - 
Dunne polynomials of genus one, one should take care that the trial function (j) belongs to a 
finite-dimensional invariant subspace Wm- This can always be done for quasi-exactly solvable 
models of quantum mechanics. Indeed, consider a quantum model whose hamiltonian can be 
represented in the form 

H = Y,C^pJo^{h)J^{h) + Y,CMh) (9.1) 

where Cap and Ca are some numbers and Ja{h) are generators of a certain Lie algebra 
realising a representation with highest weight h and having the form of first-order differential 
operators. In this case, it is natural to identify the vector cj) with the highest weight vector \h). 
If the representation is infinite-dimensional and irreducible, then all the terms of the Lanczos 
sequence 0„ = Pn{H)\h) are linearly independent and this sequence is infinite. Otherwise, 
the operator H would have a finite-dimensional invariant subspace, which would contradict 
to the condition of irreducibility of the representation. Assume now that the representation 



is finite-dimensional, so that the model (|9.1j ) is quasi-exactly solvable ||5[. Since the highest 
weight vector \h) belongs, by definition, to the representation space, the Lanczos sequence 
will contain in this case only a finite number of linearly independent terms. According to 
general reasonings given above this naturally leads us to Bender - Dunne polynomials of 
genus one. 
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